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Abstract. The first determination of the age of the Galactic thin disk from Th/Eu nucleocosmochronology was accomplished by 
us in Papers I and II. The present work aimed at reducing the age uncertainty by expanding the stellar sample with the inclusion 
of seven new objects - an increase by 37%. A set of [Th/Eu] abundance ratios was determined from spectral synthesis and 
merged with the results from Paper I. Abundances for the new, extended sample were analyzed with the aid of a Galactic disk 
chemical evolution (GDCE) model developed by us is Paper II. The result was averaged with an estimate obtained in Paper II 
from a conjunction of literature data and our GDCE model, providing our final, adopted disk age Tq = (8.8 ± 1.7) Gyr with a 
reduction of 0.1 Gyr (6%) in the uncertainty. This value is compatible with the most up-to-date white dwarf age determinations 
(< 10 Gyr). Considering that the halo is currently presumed to be (13.5 ± 0.7) Gyr old, our result prompts groups developing 
Galactic formation models to include an hiatus of (4.7 ± 1.8) Gyr between the formation of halo and disk. 
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1. Introduction 

The age of the Galactic thin disk 1 is an important constraint for 
Galactic formation models. It is usually estimated by dating the 
oldest open clusters with isochrones or white dwarfs with cool- 
ing sequences. These methods are strongly dependent on stellar 
evolution models and on numerous physical parameters known 
at different levels of uncertainty. Nucleocosmochronology is 
only weakly dependent on main sequence stellar evolution 
models, allowing for a nearly independent crosscheck of other 
techniques. 

We were the first to determine an a ge for the Galactic disk 
from Th/Eu nucleocosmochronology - del Peloso et al 1l2005bl 



Send offprint requests to: E.E del Peloso 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile, under the ESO programs and the ESO- 
Observatorio Nacional, Brazil, agreement, and at the Observatorio 
do Pico dos Dias, operated by the Laboratorio Nacional de 
Astrofi'sica/MCT, Brazil. 

** Full Table [5] is only available in electronic form at the CDS via 
anonymousftptocdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? 

1 All references to the Galactic disk must be regarded, in this work, 
as references to the thin disk, unless otherwise specified. 



(Paper I) and Idel Peloso et alJl2005al (Pap er II), based on the 
PhD thesis of one of us ( Idel Pelosdl20 03). This work, the last 
part in a series of three articles, aims at reducing the uncertainty 
in this determination by expanding the stellar sample from 
Pape rs I and II with objects observed only after the publication 
of del Peloso ( 2003) and deriving a new age from this extended 
sample. With this intent, we determined [Th/Eu] abundance ra- 
tios for a sample of Galactic disk stars and employed Galactic 
disk chemical evolution models developed by us in the chrono- 
logical analysis. 2 

2. Sample selection, observations and data 
reduction 

Sample selection, observations and data reduction were carried 
out following exactly the same procedures described in detail 
in Paper I. In what follows we provide a brief overview of these 
topics. 

The stellar sample of this work is composed of seven F8- 
G5 dwarfs and subgiants (Table As the objective of this 

2 In this paper we obey the customary spectroscopic notation: abun- 
dance ratio [A/B] = log 10 (Ar A //Y B ) stal - log 10 (/V A //V B )o, where JV A 
and are the abundances of elements A and B, respectively, in 
atoms cirr 3 . 
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Table 2. Radial velocities (RV) and spatial velocity compo- 
nents (U, V, and W) of the sample stars, in a right-handed 
Galactic system and relative to the LSR. All values are in 
km s _1 . 



HD 


RV 


U 


V 




W 


1461 


-9.3 


-15.3 +0.7 


-36.7 


±1.0 


+ 17.0 


±0.4 


157 089 


-161.1 


-156.0+1.1 


-35.5 


±1.0 


-2.8 


±1.8 


162 396 


-16.4 


-13.3 ±0.4 


-5.5 


±0.7 


-24.0 


±1.3 


189 567 


-10.0 


-59.7 ±1.6 


-26.4 


±0.7 


-42.6 


±2.0 


193 307 


+ 18.6 


+47.4 ±1.2 


-43.3 


±1.6 


+40.8 


±2.0 


196755 


-52.7 


-48.0 ±0.8 


-29.5 


±0.6 


-11.2 


±1.0 


210918 


-18.6 


-37.4 ±0.9 


-86.9 


±1.6 


-1.8 


±0.7 
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Fig. 1. V vs. [Fe/H] diagram for the sample stars. The diagonal 
line is a cut off between the halo and the disk populations (be- 
low and above, respectively). Velocity components V were cal- 
culated by us. The metallicities initially employed, represented 
as open squares, were taken from the literature (Table|5J>i small 
filled circles represent metallicities derived in this work. 



work is the determination of the age of the Galactic disk, we 
performed a kinematic analysis of the objects to help ensure 
that they do not belong to the Galactic halo. We calculated the 
objects' U, V, and W spatial velocity components (Table[5} and 
constructed a V vs. [Fe/H] diagram (Figure Q using metallic- 
ities fr om the literature (Table [3}. According to ISchuster et alJ 
l l 1993b . objects located above the displayed cut-off line belong 
to the Galactic halo, whereas those located below the line be- 
long to the halo. It can be seen that all sample stars are located 
far above the cut-offline. For a star to cross the line, the litera- 
ture metallicities would have to have been overestimated by at 
least 1 .2 dex, which is very unlikely. After having determined 
our own metallicities (Table|4}, we confirmed that results from 
the literature agree with our values to 0. 1 dex. 

High resolution, high signal-to-noise ratio spectra were 
obtained for all objects with the Fib er-fed Exten ded Range 
Optical Spectrograph (FEROS; iKauferet al.1l 19991) fed by the 
1.52 m European Southern Observatory (ESO) telescope, in 
the ESO-Observatdrio Nacional, Brazil, agreement (March and 
August 2001). Spectra were also obtained with a coude spectro- 



Table 4. Adopted atmospheric parameters, including the pho- 
tometric and Ho- effective temperatures used to obtain the 
adopted mean values, and the stellar masses used to obtain the 
surface gravities. 



HD 




T e „ (K) 


m 




[Fe/H] 


£(km s _1 ) 




Phot. 


Hey 


MEAN 










1461 


5727 


5705 


5717 


1.01 


4.33 


+0.17 


1.20 


157 089 


5827 


5742 


5785 


0.93 


4.09 


-0.47 


0.90 


162 396 


6072 


5979 


6026 


1.04 


4.08 


-0.37 


1.36 


189 567 


5704 


5715 


5710 


0.90 


4.36 


-0.22 


1.01 


193 307 


5999 


5956 


5978 


1.05 


4.11 


-0.24 


1.25 


196755 


5665 


5613 


5639 


1.58 


3.70 


+0.04 


1.42 


210918 


5733 


5708 


5721 


0.95 


4.27 


-0.09 


1.15 



graph fed by the 1.60 m telescope of the Observatorio do Pico 
dos Dias (OPD), LNA/MCT, Brazil (May and October 2000; 
May, August, and October 2002), and with the Coude Echelle 
Spectrometer (CES) fiber-fed by ESO's 3.60 m telescope 
(August 2003). FEROS spectra are reduced automatically 
during observation by a script executed using the European 
Southern Observatory Munich Image Data Analysis System 
(ESO-MIDAS) immediately after the CCD read out. OPD and 
CES spectra were reduced by us using the Image Reduction and 
Analysis Facility (IRAF 3 ), following the usual steps of bias, 
scattered light and flat field corrections, and extraction. 

3. Atmospheric parameters 

A set of homogeneous, self-consistent atmospheric parameters 
was determined for the sample stars, following the procedure 
described in detail in Paper I. Effective temperatures were de- 
rived from photometric calibrations (Table [3} and Ha profile 
fitting; surface gravities were obtained from T e ff , stellar masses 
and luminosities; microturbulence velocities and metallicities 
were obtained from detailed, differential spectroscopic analy- 
sis, relative to the Sun, using EWs of Fe i and Fe n lines. The fi- 
nal, adopted values of all atmospheric parameters are presented 
in Table E| 

4. Abundances of contaminating elements 

Eu and Th abundances were determined by spectral synthesis 
of the Eu n line at 4129.72 A and of the Th n line at 4019. 13 A, 
respectively. In the synthesis calculations, abundances of the el- 
ements that contaminate these spectral regions (Ti, V, Cr, Mn, 
Co, Ni, Ce, Nd, and Sm) were kept fixed in the values deter- 
mined using EWs and the atmospheric parameters obtained by 
us; see Paper I for a full description of the employed method. 
Table[5]presents a sample of the EW data. Its complete content, 
composed of the EWs of all measured lines, for the Sun and all 
sample stars, is only available in electronic form at the CDS. 
Column 1 lists the central wavelength (in angstroms), Col. 2 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Selected stellar sample. 
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I 1 1 \ 


IlK 


ri i r in ame 


K. A. 




Parallax 


V 


Spectral type 








zuuu.U 


ZUUU.U 


(mas) 




and 








(h m s) 


(d m s) 






luminosity class 


1461 


72 


1499 


00 1 8 42 


-08 03 1 1 


42.67 


6.46 


G5 V 


157 089 




84905 - 


17 21 07 


+01 26 35 


25.88 


6.95 


GO-2 V 


162 396 


6649 


87 523 - 


17 52 53 


-41 59 48 


30.55 


6.20 


F8 IV- V 


189 567 


7644 


98 959 - 


20 05 33 


-67 19 15 


56.45 


6.07 


G3 V 


193 307 


7766 


100412 - 


20 21 41 


-49 59 58 


30.84 


6.27 


GOV 


196755 


7896 


101916 k Del 


20 39 08 


+ 10 05 10 


33.27 


5.05 


G5IV 


210918 


8477 


109 821 - 


22 14 39 


-41 22 54 


45.19 


6.23 


G5 V 



References: Coordinates: SIMBAD (FK5 system); parallaxes: Hipparcos catalogue ( Perrvman & ESA 1997); visual magnitudes: Bright Star 
Catalogue l Warren & Hoffleit 1987) f or stars with an HR number an d SIMBAD for those without; spectral types and luminosity classes: 
Michigan Catalogue of HD Stars faouk & Cowle^fl975t lHouklll97i: iHouk & Swiftlll999h for all stars, with the exception of HD 196755 
(Bright Star Catalogue). 

Table 3. Photometric indices and metallicities for all sample stars. Metallicities were taken from the literature for use in the 
kinematic characterization of the sample and as input for the first step in the iterative determination of atmospheric parameters. 



HD 


(fl-V) 


Ref. 


(b-y) 


Ref. 


P 


Ref. 


(V-K) 


Ref. 


(S T - Vt) 


Ref. 


[Fe/H] 


Ref. 


1461 


0.68 


1 


0.422 


3 


2.596 


4 


1.518 


7 


0.764 


9 


+0.13 


10 


157 089 


0.57 


2 


0.379 


3 


2.584 


5 


1.434 


7 


0.619 


9 


-0.54 


11 


162 396 


0.54 


1 


0.347 


4 


2.612 


4 






0.563 


9 


-0.37 


10 


189 567 


0.64 


1 


0.406 


3 


2.583 


6 


1.516 


7 


0.718 


9 


-0.17 


10 


193 307 


0.55 


1 


0.365 


3 


2.604 


5 






0.617 


9 


-0.34 


10 


196755 


0.72 


1 


0.432 


3 










0.765 


9 


+0.01 


10 


210918 


0.65 


1 


0.404 


3 


2.590 


6 


1.500 


8 




9 


-0.12 


10 



Ref erences: 1 -Iw arren & Hoffleit! 1987); 2 - Mermilliodl 1987); 3 - Olse r] <1994 : 4 -Ejr0nbech & Olsenl 1976l ll977h: 5 -loTsed h983h : 
6 - ISchuster & NisserJ Il988h : 7 - Tdi Benedetto! Il998l) : jT HKoornneed ll983h : 9 - IPerrvman & ESAl <1997F ~10 - iTavloJ <2003l) : 11 - 
ICavrel de Strobel et alj COOlh . 



Table 5. A sample of the EW data. The complete content of 
this table is only available in electronic form at the CDS. For a 
description of the columns, see text (Sect.QJ. 



A (A) 


Element 


*(eV) 


logg/ 


Sun • 


• HD 210918 


5668.362 


Vi 


1.08 


-0.920 


8.5 • 


11.9 


5670.851 


Vi 


1.08 


-0.452 


21.6 ■ 


23.7 


5727.661 


Vi 


1.08 


-0.657 


9.5 • 


12.1 


5427.826 


Fen 


6.72 


-1.371 


6.4 • 


0.0 


6149.249 


Fen 


3.89 


-2.711 


40.9 • 


39.5 



gives the element symbol and degree of ionization, Col. 3 gives 
the excitation potential of the lower level of the electronic tran- 
sition (in e V), Col. 4 presents the solar log gf derived by us, 
and the subsequent columns present the EWs, in mA, for the 
Sun and the other stars, from HD 1461 to HD 210 918. 

Abundance results are presented in Table[6] No detailed un- 
certainty assessment was carried out. We have rather adopted 
an average of the uncertainties presented in Paper I: 0.08 dex 
for Mn, 0.09 dex for Fe, Ti, and Co, 0.10 dex for V, Cr, Ce, 
and Nd, and 0.11 dex for Ni and Sm. These values are used as 
error bars in Fig.|2j which shows the abundance patterns for all 
elements. Note that the abundances for the sample of this work 
(filled squares) agree very well with those from Paper I (open 
squares). 



Table 7. [Eu/H], [Th/H], and [Th/Eu] abundance ratios. 



HD 


[Eu/H] 


[Th/H] 


[Th/Eu] 


1461 


+0.09 


+0.08 


-0.01 


157 089 


-0.12 


-0.28 


-0.16 


162 396 


-0.23 


-0.12 


+0.11 


189 567 


+0.08 


+0.02 


-0.06 


193 307 


-0.13 


-0.01 


+0.12 


196755 


+0.02 


+0.09 


+0.07 


210918 


+0.03 


-0.01 


-0.04 



5. Eu and Th abundances 

Our determination of Eu abundances by spectral synthesis used 
FEROS spectra. As the abundances used for age determina- 
tion in Paper I were obtained with the CES coupled to the 
Coude Echelle Spectrograph (CAT), we converted our results 
to the CAT system, using the equation derived in Sect. 5.1.1 
of Paper I (see their Fig. 15): [Eu/H] CAT+CES = -0.00462 + 
0.96480 [Eu/H] FER0S . TableQpresents the [Eu/H], [Th/H], and 
[Th/Eu] abundance ratios for our sample. As uncertainty, we 
adopted an average of the values presented in Paper I: 0.04 dex 
for [Eu/H], 0.1 1 dex for [Th/H], and 0.08 dex for [Th/Eu]. 

In Fig. El we compare o ur [Eu/H] abundance ratios to 
those from Woolf etaQ ( 1 19951 WTL95), iKoch & Edvardssonl 
(120021 KE02), and Paper I. These are the best works avail- 
able with Eu abundances for Galactic disk stars, in terms of 
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Table 6. Fe, Ti, V, Cr, Mn, Co, Ni, Ce, 
used for each abundance determination. 



Nd, and Sm abundances, relative to H. N is the number of absorption lines effectively 



HD 


[Fei/H] 


N 


[Fe n/H] 


N 


[Fe/H] 


N 


[Ti/H] 


N 


[V/H] 


N 


[Cr/H] 


N 


[Mn/H] 


N 


1461 


+0.16 


52 


+0.19 




8 


+0.17 


60 


+0.16 


28 


+0.16 


5 


+0.18 


19 


+0.21 


6 


157 089 


-0.47 


31 


-0.47 




7 


-0.47 


38 


-0.21 


20 


-0.21 


3 


-0.42 


18 


-0.66 


4 


162 396 


-0.38 


41 


-0.33 




8 


-0.37 


49 


-0.29 


18 


-0.17 


2 


-0.35 


16 


-0.41 


5 


189 567 


-0.23 


44 


-0.20 




7 


-0.22 


51 


-0.12 


22 


-0.10 


3 


-0.21 


19 


-0.29 


6 


193 307 


-0.24 


45 


-0.22 




9 


-0.24 


54 


-0.21 


19 


-0.12 


4 


-0.21 


19 


-0.33 


6 


196755 


+0.04 


45 


+0.06 




9 


+0.04 


54 


+0.03 


28 


-0.01 


5 


+0.06 


20 


-0.01 


6 


210918 


-0.09 


47 


-0.09 




7 


-0.09 


54 


-0.01 


28 


-0.04 


5 


-0.08 


18 


-0.14 


6 


































HD 


[Co/H] 


N 


[Ni/H] 


N 




[Ce/H] 


N [Nd/H] 


N 


[Sm/H] 


N 










1461 


+0.21 


8 


+0.19 


10 




+0.08 


5 


+0.13 


2 


+0.17 












157 089 


-0.33 


7 


-0.45 


7 




-0.47 


5 


-0.21 


2 


-0.42 












162 396 


-0.22 


3 


-0.32 


7 




-0.36 


4 


-0.12 


2 


-0.44 












189 567 


-0.15 


8 


-0.26 


9 




-0.17 


4 


+0.04 


2 


-0.10 












193 307 


-0.17 


6 


-0.29 


8 




-0.29 


5 


-0.08 


2 


-0.39 












196755 


+0.00 


8 


+0.03 


10 




+0.03 


5 


+0.08 


2 


-0.09 












210918 


-0.05 


8 


-0.11 


9 




-0.07 


5 


+0.07 


2 


-0.03 
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Fig. 3. [Eu/H] vs. [Fe/H] diagram for our sample stars, and pig 4 [Th/H] yf , [Fe/H] diagram for Qur sample stafs> and 
those from Paper I, WTL95, and KE02. Average error bars for thofje from paper j and MKB92. Our results are in good accord 



the three data sets are provided in the lower right corner. 



with MKB92, but show a lower dispersion. 



care of analysis and sample size. Although two of our objects 
(HD 1461 and HD 157089) exhibit considerable discrepancy 
with Paper I, our results present behavior and dispersion simi- 
lar to WTL95/KE02. 

Th spectral synthesis employed CES spectra fed by the 
3.60 m ESO telescope, whereas the abundances in Paper I 
were obtained with the CES fed by the CAT. Our results were 
converted to the CAT system, using the equation derived in 
Sects. 5.2.1 of Paper I (see their Fig. 20): [Th/H] CAT+CES = 
-0.01500 + 0.86000 [Th/H] CAT+360ra . 

In Fig.Bl we compar e our [Eu/H] abundance ratios to those 
from iMorell et alJ d 19921 MKB92) and Paper I. These are 
the best works available with Th abundances for Galactic disk 
stars, in terms of care of analysis and sample size. Our results 
are in good accord with them. 



6. Chronological analysis 

In order to estimate the age of the Galactic disk, we com- 
pared the stellar [Th/Eu] abundance ratios with curves calcu- 
lated using a Galactic disk chemical evol ution (GDCE) model. 
In this model, developed by us based onlPa gel & Tautvaisienel 
(I1995I) . it was assumed that the so-called "universality of the 
r-process abundances" is valid for second and third r-process 
peaks and can be extended to the actinides. Such extension may 
not be legitim ate, as tw o ul tra-metal-poor stars - CS 31 082- 
001 (Cavreletal. 2001 and iHill etalll2002l) and CS 30306- 



( Cavrel et al. 


2001 


( Honda et al. 


2004 



132 

Th/Eu abundance ratios much higher than expected for their 
age. This could indicate that they have been formed from mat- 
ter enriched in actinides, relatively to second r-process peak 
elements. However, it is not yet clear if CS 31082-001 and 
CS 30306-132 are merely chemically peculiar objects or if 
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Fig. 5. [Th/Eu] vs. [Fe/H] diagram for our sample stars and 
those from Paper I. Curves were calculated using our GDCE 
model for the four different Galactic disk ages shown in the 
legend. 



their discrepancies could be present in other yet unobserved 
stars. A detailed description of the GDCE model was presented 
in Paper II. 

The abundances for our sample were determined in exactly 
the same way as those from Paper I, as the objective of this 
work is to expand the sample used for chronological analysis. 
Accordingly, we merged our abundance data with those from 
Paper I, resulting in a set of 28 objects. The theoretical model 
curves and the observed abundance data are presented in Fig. [5] 

The age of the curve that best fits the observed data was 
computed by minimizing the total deviation between the curves 
and the data. 26 of the 28 objects in the merged sample were 
actually used in the determination; two objects from Paper I 
were disregarded because they are too metal-rich, falling out 
of the interval where the curves are defined. Considering that 
in Paper II 19 objects were actually employed in the analysis, 
this work accomplishes a 37% increase in sample size. An un- 
certainty for the age, related solely to the uncertainties in the 
abundances, was computed through a Monte Carlo simulation. 
It must be noted that this is only an assessment of the internal 
uncertainty of the analysis, and does not take into consideration 
the uncertainties of the GDCE model itself, which are very dif- 
ficult to estimate. The uncertainty related to the model could 
very well be the main source of age uncertainty. These proce- 
dures are fully described in Paper II. The final value obtained 
using the merged abundance data set is (8.8 + 1.8) Gyr. This re- 
sult, 0.6 Gyr larger and with an uncertainty 0.1 Gyr lower than 
that of Paper II, matches very well the estimate obtained from 
literature data and the GDCE model (8.7^ Gyr). These esti- 
mates were combined using the maximum likelihood method, 
assuming that each one follows a Gaussian probability distribu- 
tion, which results in a weighted average using the reciprocal of 
the square uncertainties as weights. The final, adopted Galactic 
disk age is 

T G = (8.8 ± 1.7) Gyr. 



7. Conclusions 

We determined [Th/Eu] abundance ratios for a sample of seven 
Galactic disk F8-G5 dwarfs and subgiants. The analysis was 
carried out in exactly the same way as that of Paper I, so that 
we could merge the data, resulting in a totally homogeneous 
extended sample of 28 objects; 26 of these were actually used 
in the nucleocosmochronological analysis. A GDCE model, 
developed in Paper II, was used in conjunction with the stel- 
lar abundance data to compute an age for the Galactic disk: 
(8.8 ± 1.8) Gyr. In Paper II, Th/Eu production and solar abun- 
dance ratio data taken from the literature were analyzed with 
our GDCE model, yielding 8.7^ ? Gyr. These two results were 
combined using the maximum likelihood method, resulting in 
FINAL Tc = (8.8 + 1.7) Gyr. 

The inclusion of seven more stars in the abundance data 
base had two main consequences: it increased the age obtained 
from our stellar data by 0.6 Gyr, rendering it more compat- 
ible with the age determined from literature data, and it de- 
creased the uncertainty of the final, adopted age by 0. 1 Gyr, i.e., 
a 6% reduction. Our result remains compatible with the most 
up-to-date white dwarf ages derived from cooling sequence cal- 
culations, which indicate a low age (< 1 Gyr) for the disk 
JOswalt et alJll995l iBergeron et alJ^997t iLeggett et alJll998t 
iKnox et alll999HHansen et alJ2002l) . Considering that the age 
of the oldest halo glob ular clusters are currently estimated at 
(13.5 ± 0.71 Gvr JPont et alJl998H.Timenezll999HGratton et alJ 
200i lKrauss & Chaboverll2003h . an hiatus of (4.7 + 1.8) Gyr 
between the formation of halo and disk must be taken into con- 
sideration in future Galactic formation models. 
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